Turbulence and Little Rip Cosmology 
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A variety of conditions is considered under which the cosmic dark fluid may be able to develop 
a future Big Rip or Little Rip singularity. Both one-component and two-component models are 
considered. In the last-mentioned case we present a way in which the fluid can be decomposed into 
(3JQf two components, one non-turbulent (ideal) and one turbulent part, obeying two different equations of 

^ , state. For the non-turbulent part, the thermodynamical parameter, commonly called w, is assumed 

to be less than —1 throughout. For the turbulent part, it turns out that it is sufficient that uwb lies 
in the quintessence region in order to lead to a singularity. Both Big Rip and Little Rip behaviour for 
dark energy are found. In the one-component case, we examine how the universe may develop from 
<n: a viscous era with constant bulk viscosity into a turbulent era, the turbulence in effect protecting the 

universe from encountering the singularity at all. The equivalent description of the same cosmology 
' in terms of inhomogeneous (imperfect) fluid is also presented. 

CT 
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6J[) I. INTRODUCTION 

It lias become customary to explain the observed acceleration of the universe [H, Hi in terms of dark energy fluid 
" (for recent reviews, see 0, I3])j which is expected to have the strange properties like negative pressure and/or negative 
, entropy. According to the latest supernovae observations the dark energy amounts to about 73% of the total mass 
t"**- ■ energy of the universe Q. Although astrophysical observations favor the standard ACDM cosmology, the equation- 
of-state (EoS) parameter w is still determined with uncertainty: it is not clear if w is less than —1, equal to —1, or 
■ larger than — 1. Current observations suggest that w = — 1.04+oao IEQ- 

The very interesting but least theoretically understood case corresponds to w < — 1 (phantom dark energy) where 
. all four energy conditions are violated. Although the theory is unstable from a quantum field theoretical viewpoint, 
it could be stable in classical cosmology. There are observations [1] indicating that the crossing of the cosmological 
constant/phantom divide took place in the near past (or will occur in the near future). An essential property of (most) 
phantom dark energy models is the Big Rip future singularity @ (see also [13; El): where the scale factor becomes 
infinite at a finite time in the future. A softer future singularity caused by phantom or quintessence dark energy 
is the sudden (Type II) singularity [l2j where the scale factor is finite at Rip time (for classification of finite-time 
singularities, see [13( ). Recently, an attempt to resolve the finite-time future singularities has been proposed in the face 
of mild phantom models where w asymptotically tends to —1 and where the energy-density increases with the time 
or remains constant, but where singularity occurs in infinite future [T3 - [l7| . The key point here is that if w approaches 
— I sufficiently fast, then it is possible to have a model in which the time required for the occurrence of a singularity 
is infinite, i.e., the singularity effectively does not happen. However, if the energy density grows, the disintegration of 
bound structures necessarily has to take place, in a way similar to the case of a Big Rip. Such Rip phenomena turns 
out to be common for Big Rip, Little Rip or Pseudo-Rip cosmologies, destroying all bound structures in a finite time. 
It is remarkable that mild phantom scenarios like Little Rip or Pseudo-Rip may easily mimic our current ACDM era 
and indicate to quite long existence (billions of years!) of our universe before the Rip occurs. 

The pioneering works on the future singularity Q were considering the cosmic fluid to be non- viscous. This is an 
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idealized model, of course; it is useful in practice in many cases but generally not able to cope with intricate micro- 
scale phenomena that occur especially near solid boundaries. A next step in complexity is to allow for deviations 
from thermal equilibrium to the first order. That means, one has to introduce two viscosity coefficients, namely the 
shear viscosity 77 and the bulk viscosity £. 

In the present paper we consider dark energy era with future Big Rip or Little Rip singularity under a variety of 
conditions, assuming a more realistic form for the cosmic fluid than what has commonly been the case. We start 
in the next section by reviewing essentials of the theory of a dark fluid satisfying the condition p < —p. i.e. a 
phantom fluid. While in the general case one has to account for for both coefficients r\ and £, as mentioned, we 
shall assume in conformity with usual practice that spatial anisotropies are smoothed out. Thus only ( has to be 
included. Viscous Little Rip cosmology in an isotropic fluid has recently been worked out; cf. Eq. ([7]) below. In Sec. 
Ill we focus attention on the description of a turbulent state of the fluid in the late universe. Such a possibility would 
seem physically most natural, in view of the violent motions expected in the vicinity of the singularity. We propose a 
two-component model in which, for physical reasons, the turbulent energy component is put proportional to the scalar 
expansion. Wc consider in this context two different proposals for the thermodynamic parameter wturb occurring in 
the equation of state. First, we put wturb equal to the usual parameter w for non-turbulent matter, meaning that the 
turbulent matter component behaves as a passive ingredient as far as the equation of state is concerned. Second, we 
allow for Wturb > — 1 5 meaning that we cover the region — 1 < itfturb < also. Both Big Rip and Little Rip evolutions 
for dark energy are found. In Sec. VI we consider an approach which is quite different, namely, we model the universe 
as a one-component dark fluid that becomes suddenly transformed into a turbulent state. A noteworthy property of 
this model is that the turbulence in effect protects the universe from encountering the future singularity at all. In 
Sec. VII we present an equivalent description in terms of an inhomogeneous imperfect fluid, including the effect of 
viscosity but not that of turbulence. Some summary and outlook is given in Discussion section. 



II. DARK FLUID WITH A BULK VISCOSITY 



Consider the following model for the cosmic fluid in the later stages of the development of the universe, assuming 
that the fluid is dark and satisfies the inequality 

P<-P, (1) 

so that the equation 

H = - 1 -k\ P + P ) (2) 
with H — a/a and k 2 = 8irG implies the property 

H>0. (3) 

This is the phantom region, corresponding to the thermodynamic parameter w = p/ p being less than —1. In the usual 
case of this type universe the future singularity is a true mathematical singularity, reached in a finite time, and is 
called the Big Rip. It is notable, however, that in the limiting case when w — > — 1 from below, the future singularity 
is only asymptotically reached. This scenario is called the Little Rip. 

It is noteworthy that in a cosmological context the value of 77 appears to be very much larger than that of £. At 
least that is so in the earlier stages of the development of the universe where the physical conditions are better known 
than in the later stages and conventional kinetic theory can be used to calculate the viscosities. A calculation of the 
viscosity coefficients was made by Caderni and Fabbri For instance, considering the instant t = 1000 s after the 
Big Bang, it turns out that 

j] = 2.8 x 10 14 g curV 1 , C = 7.0 x 10~ 3 g curV 1 , (4) 

(cf. also Ref. [l9|). showing the large difference in magnitude between r\ and £. Yet, it appears that spatial anisotropies 
in the universe are effectively smoothed out, at least on a large scale, so that in most current models the universe is 
assumed to be spatially isotropic. It means that £ is retained, but r\ omitted in the Friedmann equations. 

A theory of viscous Little Rip cosmology was recently given in Ref. HTJ] . Let us recapitulate one of the characteristic 
results from that investigation: If the effective pressure p G g is assumed to have the explicit form 



Pes = -p~ Ayfp - 3(H 



(5) 
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with A a positive constant, and if moreover the bulk viscosity £ is assumed to satisfy the condition 

3(H = £o = constant , (6) 
then the following expression is found for the time dependent energy density 



P(t) = 
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(7) 



Here subscript zero refers to the present time. It is thus an infinite time needed to reach the infinite energy density 
case. This is precisely the characteristic property of the Little Rip phenomenon. 

III. THE TURBULENT APPROACH 

Let us now apply a physical point of view on the dark energy universe, in its later stages when it approaches 
the future singularity. The simple description above, in terms of macroscopic bulk viscosity in the fluid, cannot be 
considered to be satisfactory, due to the following reason: In the assumed states of violent local fluid element motion 
a transition into turbulent motion seems to be inevitable. The local Reynolds number must be expected to be very 
high. That brings, in fact, the shear back into the analysis, not in a macroscopic sense as before, but in a local sense 
causing the distribution of local eddies over the wave number spectrum. What kind of turbulence should we expect? 
The natural choice is that of isotropic turbulence, which is a topic reasonably well understood. Thus, we should 
expect a Loitziankii region for low wave numbers where the energy density varies proportionally to k ; for higher k 
we should expect an inertial subrange characterized by the formula 

E{k) = ae 2 / 3 k- 5 / 3 (8) 

with a the Kolmogorov constant and e the mean energy dissipation per unit time and unit mass; and finally when 
the values of k become as high as the inverse Kolmogorov length t]l, 

1 / e \V4 

k^k L = — =[ — ) (9) 



with v the kinematic viscosity, we enter the dissipative region where the local Reynolds number is of order unity and 
heat dissipation occurs. In accordance with common usage we shall consider the fluid system as quasi-stationary, and 
omit the production of heat energy. In practical cases it may be useful to combine these elements into the useful von 
Karman interpolation formula which covers the whole wave number spectrum (cf., for instance, Refs. (2l| - [23l |). 

However, the full spectral theory of isotropic turbulence will no be needed in our first approach to the problem. 
Rather, we shall in the following focus attention on how the turbulent part of the energy density, called pturb, can be 
estimated to vary from present time to onwards. First, we write the effective energy density as a sum of two terms, 

PcS = P + Pturb , (10) 

where p denotes the conventional macroscopic energy density in the local rest inertial system of the fluid. It is natural 
to assume that pturb is proportional to p itself. Further, we shall assume that ptm-b is proportional to the scalar 
expansion 6 = U^-^ = 3H. This because physically speaking the transition to turbulence is expected to be more 
pronounced in the violent later stages, and a proportionality to the scalar expansion is mathematically the most 
simple way in which to represent the effect. Calling the proportionality factor r, we can thus write the effective 
energy density as 

PeS =p(l + 3rH). (11) 
Consider next the effective pressure p e ff- We split it into two terms, 

PcS = V + Pturb , (12) 

analogously as above. For the conventional non-turbulent quantities p and p we assume the standard relationship 



p = wp, 



(13) 
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where —1 < w < —1/3 in the quintessence region and w < —1 in the phantom region. The question now is: How 
does pturb depend on pturb? There seems to be no definite physical guidance to that problem, so we shall make the 
simplest possible choice in the following, namely write 

Pturb = Wt ur b Pturb , (14) 

with u>turb a constant. 

We shall consider two different possibilities for the value of wturb- The first is to put wturb equal to w in Eq. (fl"3|) . 
meaning that the turbulent matter behaves in the same way as the non-turbulent matter as far as the equation of state 
is concerned. This option is straightforward and natural, and is not quite trivial since pturb and p behave differently, 
in view of Eq. Our second option will be to assume that lUturb takes another, prescribed value. In view of the 

expected violent conditions near the future singularity, it might even be natural here to chose the value u>t U rb = +L 
i.e., the Zel'dovich fluid option. 

The first and the second Friedmann equations can now be written 

H 2 = ^k 2 P (1 + 3tH) , (15) 



— + H 2 = -n 2 p(w + 3riiu;turb) • (16) 
a 

(recall that k 2 = 8irG) . This may be compared with earlier attempt to introduce the turbulence in dark energy [24| . 

Equations (TT5|) and (TT6"]) determine our physical model. Recall that its input parameters are {w, w tU rb, t}, all 
assumed constant. From these equations we can now describe the development of the Hubble parameter. For 
convenience we introduce the quantities 7 and 7t U rb, defined as 

J=l+W, 7turb = 1 + Wturb ■ (17) 

We can then write the governing equation for H as 

(I + Zt^H+^^ + ^t^H 3 = 0. (18) 
This equation is in principle to be integrated from present time t = to = onwards, with initial value H = Hq = do/ao- 



A. On the energy balance equation 

If T t Qj denotes the total energy-momentum tensor for the cosmic fluid, we must have 

r&.„ = o, (19) 

as a consequence of Einstein's equation. 

In most cases studied, the expression for T^ t can be written down explicitly; this is so for non-viscous fluids as 
well as with macroscopic viscous fluids. In the present case this no longer true, however, since the turbulent energy 
is produced by shear stresses on a small scale, much less than the scale of the macroscopic fluid equations. That is, 
we are dealing with a non-closed physical system, of essentially the same kind as encountered in phenomcnological 
electrodynamics in a continuous medium in special relativity. It implies that the source term in the energy balance 
equation has to be put in by hand. 

Let henceforth T^ v refer to the non-viscous part of the fluid. We may express the energy balance as 

p + 3H(p + p) = -Q 1 (20) 

where the source term Q is positive, corresponding to an energy sink for the non-viscous fluid. We shall put Q equal 
to ep, where the specific energy dissipation e however shall be taken to involve the large Hubble parameter H in the 
later stage of the development. Let us assume the form 

e=eo(l + 3rff), (21) 

eo being the specific energy dissipation at present time. This equation is seen to contain the same kind of development 
as assumed before; cf. the analogous Eq. (fTTj). Thus, our ansatz for the energy balance reads 

p + 3H{p + p) = -pe (l + 3TH). (22) 
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We shall in the following consider some examples. First, however, it is of interest to compare Eq. ([22]) with the generic 
equation commonly accepted in order to deal with the extra pressure in a fluid, 

= pi^u" + ( P + n)/^" , (23) 

with W = g^ v + U^U U the projection tensor. Here IT is the extra pressure brought about by a variety of effects, 
among them typically viscosity, matter creation, or eventually a combination of these. (For a more detailed discussion 
along the lines cf., for instance, Refs. (25l - [27| .) In the case of bulk viscosity, it is known that II = —3(H. Thus, 
equation (|22p may be regarded as an energy equation analogous to the pressure equation (|23[) . 



IV. THE CASE TOturb =w<-l 



This case means that the turbulent component of the fluid is regarded as a passive ingredient as far as the equation 
of-state (EoS) parameter is concerned. The time development of p and pturb will however be different. Equation (|18| 
reduces to 



£T+^ 7 F 2 = 0, (24) 



leading to 



where we have defined 



(25) 

Z = 1 + ^H Q t (26) 

(note that w < — 1 implies 7 < 0). Thus we have a Big Rip cosmology, where the future singularity time t s is given 
by 

3|7|i/ 

The scale factor becomes correspondingly 

a = a Q Z 2 ^ , (28) 
and from the first Friedmann equation (1151) we get the non-turbulent energy density as 



3ff 2 1 1 



P k 2 Z Z + 3tH 
The ratio between turbulent and non-turbulent energy becomes 



(29) 



Pturb =3rF = ^. (30) 



p Z 

It is of main interest to consider the behavior near t s . As Z = 1 — t/t s we see that 



H ~ , a ~ (31) 

ta-t (i s -i) 2/3171 



P-j-^, (32) 

Cs — I P t s —t 

Notice the difference from conventional cosmology: the behavior of H and a near the singularity is as usual, while 
the singularity of p has become weakened. The reason for this is, of course, the non-vanishing value of the parameter 
t. Moreover, as t — > t s all the non-turbulent energy has been converted into turbulent energy. From a physical point 
of view, this is just as we would expect. 
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Let us examine how this formalism compares with our ansatz (|22j) for the energy balance of the non-turbulent fluid. 
The left hand side of that equation is explicitly calculable by means of the formulas just derived. It is convenient here 
to make use of the mathematical relationship 



by means of which obtain 



j> + 3H(j,+p)=a-*j t (pa^), (33) 



In the limit t — > t„, this reduces to 



p + 3H( P + P )^^-^M±. (35) 



In Eq. ([22]) we calculate the right hand side 



«< 1 + ^> = -^!(CT)( 1 + ^> 



which for t — > t s reduces to 



- P eo(l + 3rH)^-^-^. (37) 

This has actually the same form as Eq. (|35j) . thus supporting the physical consistency of the ansatz. We obtain the 
relationship 

€0= 2V (38) 

This could hardly have been seen in advance. The specific energy dissipation eo at the present time is related to the 
EoS parameter 7 and the parameter r introduced in Eqs. (fTT|) and (fT2"|) . in a simple way. In geometric units, the 
dimension of eo is cm -1 . 



A. A comment on Little Rip cosmology 



The theory given above concerns Big Rip, a singularity obtained in the future at a finite time. A milder variant 
of this is the Little Rip scenario, corresponding to an infinite span of time needed to reach the singularity (cf., for 
instance, Refs. [28| and [20j). To achieve a Little Rip, we have to modify to some extent the above basic assumptions. 
One such modification is the following: 

1. Take the equation of state to be 

p = - P -A^p, (39) 

with A a positive constant. 

2. Put t = in the first Friedmann equation (|T5|) . 

3. Assume a milder form for the sink term in the energy balance equation than the form given in (|22j) . A simple 
option is to take Q to be proportional to H itself. Calling the proportionality constant £0 (> 0), we get the 
following modified ansatz 

p + 3H(p + p) = -SoH. (40) 

From Eqs. (|3^|) and (|4T))) . we now have 



p = (ZAJ-p - Hv)H . (41) 
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Comparing this with the modified first Fricdmann equation we get 

V3 f dp 2 1 



t = 



Inverting this equation we have 



In 



y/p(3Ay/p-£ ) \/3kA 3Ay/n-$ 



c2 

9A 2 



1 



- 1 



exp 



(42) 



(43) 



This expression shows just the characteristic property of the Little Rip phenomenon: the universe reaches the 
state p — > oo, but needs an infinite time to do so. 

One may ask: Are the physical assumptions underlying the Big Rip scenario stronger than those underlying the Little 
Rip one? In our opinion this is most likely so, although there are of course considerable uncertainties connected with 
the future time development of the universe. 



V. THE CASE w < -1, w tUTh > -1 



This case is thermodynamically quite different from the preceding one as the turbulent component of the fluid is no 
longer a passive ingredient. We have now 7t U rb = 1 + Wturb > 0, which means that we cover the region — 1 < tfturb < 
also. In the latter region, the turbulent contribution to the pressure is still negative as above, while if Wturb > the 
turbulent pressure becomes positive, just as in ordinary hydrodynamical turbulence. 

The governing equation (|18|) . written as 



(1 + 3tH)H = -H 2 (\ 1 \ 



3T7t urb -H") 



tells us that at the present time t = the condition 



\l\ > 3T7turbH"o 



(44) 



(45) 



must hold. This is so because at the present time, the turbulent part is regarded as unimportant. This corresponds 
to the inequality H > (cf.. Eq. (J3J)). 

Equation (|44| can be integrated to give t as a function of H, 



t = 



3|7l 



1 1 \ 2r / 7tuxb 
#o H) | 7 | V | 7 | 



In 



M ~ 3T7turb# H 

| 7 | - 3r7turb-ffo H 



A striking property of this expression is that it describes a Little Rip scenario. As t 
reaches a finite critical value 



H CI it 



1 



17 1 



(46) 

oo, the Hubble parameter 
(47) 



3r 7tur b 

The physical role of 7 t U rb is thus to postpone and weaken the development towards the future singularity. 

A natural choice for the EoS parameter u>t U rb m the vicinity of the singularity, in view of the violent motions 
expected, would be 

Wturb = +1 , (48) 

that means, a Zcl'dovich fluid. This is an extreme case, where the velocity of sound equals the velocity of light. 



A. Vacuum non-turbulent fluid component (w — —1), i«turb > — 1 

This case, corresponding to a vacuum non-turbulent fluid (p = —p), has to be considered separately. From Eq. (|44[) 
we can solve directly for H as a function of t, 



H_ = StHq + y/1 + 6tH + 9rffg(T + 7turb t) 

H l + 6TH + 9TH% ltulh t ' [ ' 

Thus as t — > oo, H decreases smoothly to zero as t -1 / 2 . We conclude that at least quintessence conditions are 
necessary to produce a future singularity. 
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VI. A ONE-COMPONENT DARK FLUID 

We now turn to an approach that is quite different from the one above, namely to consider the cosmic fluid as 
a one-component fluid. Thus the distinction between a non-turbulent and a turbulent fluid component is avoided 
altogether. This new approach is actually more close to the usual picture in hydrodynamics, where a fluid is known 
to shift suddenly from a laminar to a turbulent state. 

Consider the following picture: the universe starts from present time t = as an ordinary viscous fluid with a bulk 
viscosity called £, and develops according to the Fricdmann equations. We assume as before that the EoS parameter 
w < — 1, meaning that the universe develops in the viscous era towards a future singularity. Before this happens, 
however, at some instant t = i*, wc assume that there is a sudden transition of the whole fluid into a turbulent state, 
after which the EoS parameter is wturb and the pressure accordingly pturb = Wturb Pturb- As before we assume that 
ifturb > — 1, and for simplicity we take as well as w and wturb, to be constants. One may ask: What is the resulting 
behavior of the fluid, especially at later stages? 

The problem can easily be solved, making use of the condition that the density of the fluid has to be continuous at 
t = t*. In the viscous era < t < the energy- momentum tensor of the fluid is 



where 



= pU^Uv + (p - 3#C) V . (50) 



= sv + U^U V (51) 



is the projection tensor (the shear viscosity is omitted because of spatial isotropy). Solving the Friedmann equations 
one gets H3| 

H e t/tc 

H = i 3 ,\>y\H t c (et/t<-l)' (52) 



2 

a 



[l-f|7|#o*c(e^-l)] 

' [l-§| 7 |ff t c ( e V*c_l) 



2/3|-7| 



(53) 



(54) 



where t c is the 'viscosity time' 



(^ 2 C) • (55) 



The values H*, a*, at t = i* are thereby known. 

In the turbulent era t > i* we can make use of the same expressions ([5^1) - (|5"4")) as above, only with substitutions 
t c — >• oo (( —> 0), t —> t — t*, w — > Wtmb, H — >• H*. a — >• a*, po —> p*. Thus 

H = o r , (56) 

l + §7turbff* (t-U) 

(57) 



[l + f7tu rbJ ff*(^-^)] 2/37turb ' 
P = 2 (58) 

[l + §7turb#» (*-**)] 

(recall that 7t U rb > 0). Thus the density p, at first increasing with increasing t according to Eq. (pyfj). decreases again 

once the turbulent era has been entered, and goes smoothly to zero as t~ 2 when t — > oo. In this way the transition 

to turbulence protects the universe from entering the future singularity. 

It should be noted that whereas the density is continuous at t = t* the pressure is not: In the laminar era 
= wp* < 0, while in the turbulent era = Wtm-b P* wm " even be positive, if ifturb > 0. Thus, we demonstrated the 

possible role of turbulence to protect the universe from the future singularity. In the same fashion, one can consider 

its role in protecting the universe from Rip, i.e. disintegration of bound structures. 
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VII. INHOMOGENEOUS (IMPERFECT) DARK FLUID DESCRIPTION 

Let us present the equivalent formulation in terms of inhomogeneous (imperfect) fluid. By including the effect by 
the viscosity but not including that by the turbulence, we may consider the following EoS 

P = -p + f(p)-KH. (59) 

Here f(p) is an appropriate function of the energy-density p. The bulk viscosity £ can be a function of p and H, 
C = C (p,H). Let assume the Hubble rate H is given by a function H = h(t). Then by using the first FRW equation, 
we find 

t = h-\H) = h' 1 (k*/£) . (60) 
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Then by using the second FRW equation, we obtain 

f{p) = ~H + 3^(p,H)H = ~h , (h- 1 (K^^+Ky/^(p,K^ . (61) 

Conversely, if f(p) is given by (|6"Tj) . we obtain a solution of the FRW equation as H = h(t). 
Just for the simplicity, we may assume iC,H is a constant as in 

K{p, H)H = K^TpC (p, K^j = Co • (62) 

As a simple example, we may consider the model (|25[) with (|26l) . which is realized by including the turbulence. Then 
we obtain 

/(p)=7P + Co. (63) 
As another example, we consider the model (|44]l , which give (|45j) . Since 

rr 3g 2 (| 7 | - 3r 7turb g) k 2 P (|t| - r 7turbK V37) 

h ®= H = 2(1 + 3tH) = 2(1 + ™^) ' (64) 

we find 

f(P) = 77— 7==v + Co • (65) 

(1 + TKsfZp) 

Then the models with turbulence can be equivalent to the inhomogeneous (imperfect) fluid models without turbulence 
with respect to the expansion history of the universe. In other words, the effect of turbulence and/or of viscosity may 
be always traded via the change of the effective equation of state. 



VIII. DISCUSSION 



In summary, we have studied a dark fluid universe where a finite-time future Big Rip or a Little Rip cosmology 
occurs. Special attention is paid to the role of viscosity via its elementary physical properties. The possibility of a 
viscous Little Rip cosmology where singularity occurs in the infinite future is confirmed. We propose to take into 
account a turbulent state for the dark fluid in the late universe. An explicit two-component model with turbulent 
component proportional to the Hubble rate is developed. A conventional as well as a quintessence value for the 
turbulent component equation of state parameter is considered, and the occurrence of a Big Rip and a Little Rip 
cosmology is demonstrated. However, for a one-component fluid which suddenly transforms into a turbulent state, 
the late-time acceleration universe is found to be qualitatively different: the turbulence may protect the universe from 
encountering a future singularity at all. It is interesting that an equivalent description in terms of an inhomogeneous 
fluid may also be developed. 

Turbulence is a fundamental property of classical fluids. It is known to be important not only theoretically but also 
for number of practical applications in everyday life (the air plane flight turbulence is a well known example). The 
mysterious dark energy is often considered as some kind of classical fluid with unusual properties. It is then natural 
to expect (if such a picture is correct) that turbulence phenomena may be important also for dark energy especially 
in the very late violent universe. In this work we have made a first step towards the construction of a bridge between 
turbulence and dark energy. It turns out that turbulent dark energy may have quite rich properties, in particular, in 
predicting the possible absence of future singularities. However, only forthcoming precise observational data about 
dark energy may confirm the role of turbulence in the dark energy paradigm. 
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